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On the geometry of entangled states. 
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The basic question that is addressed in this paper is finding the closest separable state for a given 
entangled state, measured with the Hilbert Schmidt distance. While this problem is in general very 
hard, we show that the following strongly related problem can be solved: find the Hilbert Schmidt 
distance of an entangled state to the set of all partially transposed states. We prove that this latter 
distance can be expressed as a function of the negative eigenvalues of the partial transpose of the 
entangled state, and show how it is related to the distance of a state to the set of positive partially 
transposed states (PPT-states). We illustrate this by calculating the closest biseparable state to 
the W-state, and give a simple and very general proof for the fact that the set of W-type states 
is not of measure zero. Next we show that all surfaces with states whose partial transposes have 
constant minimal negative eigenvalue are similar to the boundary of PPT states. We illustrate this 
with some examples on bipartite qubit states, where contours of constant negativity are plotted on 
two-dimensional intersections of the complete state space. 
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In this paper we try to get some insight into the ge- 
ometrical structure of entangled states. The main goal 
will be to characterize the distance of an entangled state 
to the set of separable states. Related questions were ad- 
dressed in the papers of Zyczkowski et al. PH3], Pittenger 
et al. Q and Witte et al. || (see also Ozawa ||), although 
here we attack the problem from a different perspective. 

The concept of negativity will turn out to be very much 
related to the Hilbert-Schmidt distance of a state to the 
set of separable states. It originates from the observa- 
tion due to Peres M that taking a partial transpose of a 
density matrix associated with a separable state is still 
a valid density matrix and thus positive (semi)dcfinitc. 
Subsequently M.Horodecki,P.Horodecki and R.Horodecki 
H proved that this was a necessary and sufficient con- 
dition for a state to be separable if the dimension of the 
Hilbert space does not exceed 6. In higher dimensional 
systems, no easy way of determining the separability of 
a state exists due to the existence of bound entangled 
states. We will therefore content ourselves to calculate 
the Hilbert Shmidt distance of an entangled state to the 
set of PPT-states (Remark that in the case of two qubits 
no bound entangled states exist). 

This problem is highly related to calculating the dis- 
tance of an entangled state to the set of partially trans- 
posed states, as the intersection of the set of all states 
with the set of all partially transposed states is equal 
to the set of all PPT-states. This is visualized in figure 
(|l|), where the boundary of the convex set of states H 
consists of rank deficicient states. The set of partially 
transposed states is completely isomorf with the set of 
states, and can be seen as some kind of reflection of the 
set of states. The intersection of both sets is the convex 
set of PPT-states. 

From figure (|l|) it is immediatly clear that the dis- 
tance of an entangled state to the PPT ones is equal to 
the distance of an entangled state to the set of partially 



transposed states iff the closest partially transposed state 
is positive (semi)-definite; this condition will turn out to 
be almost always true. 




FIG. 1. The set of all states is depicted by H and the set 
of all partial transposed states by H PT . The intersection of 
both is the set of all PPT-states. 

Let us now calculate the closest partially transposed 
state to an entangled state. The key observation is the 
fact that the Hilbert Schmidt norm is preserved under the 
partial transpose map. Therefore the proposed measure 
can be defined in the space of partial transposed density 
matrices as the minimal Hilbert Schmidt distance of p PT 
to the surface of positive (semi)-definite matrices with 
trace 1, this surface being the partial transpose of the 
boundary of PPT-states. 

We are therefore looking for the best positive semidef- 
inite approximation of the indefinite matrix p PT in the 
Hilbert-Schmidt norm: 



l|p PT -^l| 2 = v Tr((p^- f T)2) 
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Writing the eigenvalue decomposition p PT — UDU\ and 
absorbing U in a, this problem is equivalent to finding 



a such that \\D — a\\ is minimal. Using the eigenvalue 
decomposition a = VE 2 V J( with Tr (E 2 ) — 1, this can 
be written as a Lagrange constrained problem with cost 
function: 



K =\\D- VE 2 V^\\ 2 - A (Tr (E 2 ) - l) 



(2) 



It is immediatly clear that the optimal unitary V is given 
by the identity: a positive definite matrix remains posi- 
tive definite if off-diagonal elements are made zero. Dif- 
ferentiation leads to the result that the e 2 are either equal 
to 0, either equal to d, + A. Normalization fixes the 
value of A. Straightforward calculations show that th 
e 2 corresponding to the negative eigenvalues dj have to 
be choosen equal to zero and the other ones either equal 
to di + A either equal to 0, depending on the sign of dt + A. 
The algorithm for finding the closest partially transposed 
state therefore becomes: 

1. Calculate the eigenvalue decomposition of p PT = 

udW 

2. Define E 2 as the unique diagonal positive (semi)- 
definite normalized matrix such that its elements 
are ef = d,- + A or ef = 0. 

3. The closest partially transposed state p s is given by 

p s = (UE 2 W) . The Hilbert Schmidt distance 
between both states is given by 



\P-Ps\\2 = 
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where I n is the set of all indices corresponding to 
the negative eigenvalues of p PT , I p is the set of in- 
dices corresponding to positive eigenvalues of p PT 
but for which E 2 
E 2 . 



0, and rip denotes the rank of 



If p s is a state, it is guaranteed to be the closest PPT- 
state. Numerical investigations show that for example in 
the two qubit case the positiveness of p s happens in ap- 
proximately 97% of the cases. If p s is not positive, then 
the distance to the set of partially transposed states cal- 
culated is a (fairly good) lower bound on the distance of 
the entangled state to the set of PPT-states. 

Let us illustrate the above procedure with an exam- 
ple. Say we want to find the closest biseparable 2x4 
state to the three qubit W-state f(j) \W) = (|001) + 
|010) + 1 100))/v / 3- The eigenvalue decomposition of 
(|VK)(M / |) PT = UDU\ with the partial transpose opera- 
tion taken over the 4-dimensional Hilbert space, is given 
by: 



D = diag( 2/3 v^/3 1/30000 -a/2/3 ) 
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l/v/2 . . l/v/2 
1/n/2 . . -1/V2 
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The eigenvalues E 2 are readily obtained: 

E 2 = diag ( 2/3 - x/2/9 2^/9 1/3 - \/2/9 ) 

(6) 

Taking the partial transpose leads to the state p s , where 
we used the notation c = V2/18: 
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The eigenvalues of p s are non-negative and it is possi- 
ble to show that p s is separable. We have therefore 
found the closest biseparable 2x4 state to the \W)- 
state, and the Hilbert Schmidt distance to it is equal to 
(2/3) . Recently, the question arised whether the set 
of W-type states is of measure zero. Using the language 
of the Hilbert-Schmidt distance, this problem is readily 
solved. Indeed, the question is solved if we can prove 
that the state obtained by mixing the W-state with a 
small random completely separable mixed state remains 
outside the set of all convex combinations of bisepara- 
ble states (with relation to whatever partition) . As there 
is no biseparable pure state infinitesimally close to the 
W-state, and a mixed state not infinitesimally close to 
a pure state is always at a finite distance from whatever 
pure state, it is proved that the set of W-type states is 
indeed not of measure zero. A different proof was given 
by Acin et al. [Q. Remark that the above proof is very 
general and can be used in systems of arbitrary dimen- 
sions: whenever there exists a pure state ipi that can 
probabilistically be converted into another one ■02 but 
not vice-versa, the set of ^»i-like states minus the set of 
the -02-like states is of finite measure if there does not 
exist a T/^-hke state infinitesimally close to ij)-\\ 

The concept of negativity is also connected to the con- 
cept of robustness of entanglement [Q. Indeed, let us 
calculate how much an entangled bipartite state of what- 
ever dimension has to be mixed with the identity before 
it gets PPT. In analogy with the previous derivation of 
the Hilbert-Schmidt distance, this amounts to the equiv- 
alent problem of how much one has to mix the partial 
transpose of p with the identity before it gets positive 
semi-definite: 



min(l - t)p PT + -h > 
t 4 

This problem is readily solved, and the solution is 

\U>min I 



t 



(8) 
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where d m i n is the minimal negative eigenvalue of p* 
The minimal t is therefore only a function of the nega- 
tive eigenvalues. A geometrical implication of this fact 
is that all surfaces of constant d m i n are similar to the 
boundary of separable and entangled states: the set of 
all states with constant d m i n can be generated by ex- 
trapolating all lines from the identity to the boundary of 
separable states such that the distance of the extrapo- 
lated state to the identity is a constant factor (> 1) of 
the distance of the separable state to the identity. 

Let us now move to the case of two qubits. In this 
case p PT has at most one negative eigenvalue M. Nu- 
merical investigations indicate that in a vast majority of 
the states the optimal rank of E 2 is equal to three, and if 
the rank is equal to two it implies that p s has a negative 
eigenvalue. For the states for which E 2 is rank 3, it fol- 
lows that their distance to the set of partially transposed 
states is given by 

2 , . 



\\P- Ps\\ 
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where d m i n is the negative eigenvalue of p PT . Surfaces 
of two-qubit states with constant negativity, defined as 
N = 2|d m i„|, have therefore two distinct properties: they 
are all similar to each other and the Hilbert-Schmidt dis- 
tance between them is almost everywhere constant. 

Let us illustrate the above findings by explicitely cal- 
culating some two-dimensional intersections of the set 
of all bipartite qubit states including the maximally 
mixed state. In the following figures we use the met- 
ric based on the Hilbert-Schmidt distance \\pi — P2W 2 — 
Tr ((pi — P2) (pi — P2)), and directions represented or- 
thogonal to each other are orthogonal in the sense that 
Tr (A1A2) = 0. Rank deficient density operators always 
lie on the boundary of the intersection. 

Note that an explicit parameterization of the bound- 
ary between the entangled and separable states can easily 
be obtained: it is at most a quartic function of the mix- 
ing parameters of the states, as their analytic expression 
can be obtained by setting the determinant of the partial 
transpose equal to zero. 

As a first example we consider the plane containing the 
maximally mixed state and the states 
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-0.4 -0.2 0.2 0.4 0.6 0.8 
FIG. 2. Intersection of the convex set of all states includ- 
ing states (hlh and the maximally mixed state. The contours 
represent surfaces of constant negativity, the starred line is 
the boundary between separable and entangled states. 



The plane is plotted in figure (g) and the boundary of 
all (rank-deficient) states is given by the solid enveloppc. 
The starred line is the boundary between the convex set 
of separable states and the convex set of all states. The 
surfaces of constant negativity are indeed all similar to 
this boundary. The fact that the distance between these 
surfaces is not constant throughout the picture indicates 
that the closest separable states lie in other planes. Note 
that the Werner states lie along the line between the max- 
imally mixed state and the maximally entangled state p\ . 
The thirth extremal point in the undermost left corner is 
given by the rank 2 state 
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This state is called a quasi-distillable state and has some 
remarkable properties: a single copy of it can be des- 
tilled infinitesimally close to the singlet state |13n4], it 
is the state with minimal negativity for given entangle- 
ment of formation H , and it has furthermore the strange 
property that no global unitary operation can increase 
its entanglement p5| |. 

Let us now consider a different plane including the 
maximally mixed state and 
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This plane is obtained by rotating the previous plane 
around the axis pi — Ii. In this case (pi — I4) is orthog- 
onal to (p2 — Ii), and a completely different picture is 
obtained as shown in figure (£ 




The surfaces of constant negativity become straight 
lines, implying that the closest separable states lie in the 
same plane: the Hilbert-Schmidt distance betweem the 
surfaces of constant negativity has to be constant if they 
consist of parallel planes. Using the procedure previously 
outlined, it is indeed trivial to check that the separable 
state closest to the maximally entangled state p\ lies in 
the defined plane and is given by 



/ 
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(15) 



Let us rotate the plane further over the (pi — i4)-axis: 



-0.2 0.2 0.4 0.6 0.8 

FIG. 3. Intersection of the convex set of all states including 



states (fL3) and the maximally mixed state. 



Further rotation of the plane leads to the following 

states: 




" (10 10) 

mi I \ ' 



(16) 



Pi 




1 I 1 



; (0 1 1 







P2 = \\ J |(o 1 0) 



The resulting figure (||) combines the features of the 
previous figures. Three entangled disconnected regions 
arise, and once more we observe the strange shape of the 
boundary between entangled and separable states. 



(14) 



The intersection of the state space by this plane is shown 
in figure (0). 
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Intersection of the convex set of all states including 
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FIG. 4. Intersection of the convex set of all states including 



states (ll4J) and the maximally mixed state. 



FIG. 5 
states (fL6() and the maximally mixed state 



A plane with a highly symmetric contour lines is ob- 
tained if p\ and p2 are both taken to be maximally en- 
tangled states: 
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Indeed, only straight lines are obtained in figure (g). 
The thirth extremal state is in this case given by p — 
diag[l/2;0;0;l/2]. 
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FIG. 6. intersection of the convex set of all states including 
states (n7|) and the maximally mixed state. 



At last, we choose two random planes through the 
maximally mixed state and plot them in figure (pi). 





-0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 

FIG. 7. Contour plots of the negativity on random planes 
including the maximally mixed state. 



The similarity of all planes with constant negativity is 
clearly illustrated. 

We acknowledge interesting discussions with K. Au- 
denaert, A. Pittenger and K. Zyczkowski, who gave the 
idea of considering 2-D intersections of the state space. 
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